Localization of superconductivity in superconductor-electromagnet 
hybrids 



G W Ataklti 1 , A Yu Aladyshkin 1 2 , W Gillijns 1 , I M Nefedov 2 , J Van de Vondel 1 , 
A V Silhanek 1 ' 3 , M Kemmler 4 , R Kleiner 4 , D Koelle 4 , V V Moshchalkov 1 

W INPAC - Institute for Nanoscale Physics and Chemistry, Nanoscale Superconductivity and Magnetism and 
Pulsed Fields Group, K.U. Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium 

( 2 ) Institute for Physics of Microstructures, Russian Academy of Sciences, 603950, Nizhny Novgorod, 
GSP-105, Russia 

( 3 ) Departement de Physique, Universite de Liege, Bat. B5, Allee du 6 aout 17, B — 4000 Sart Tilman, Belgium 

( 4 ) Physikalisches Institut — Experimentalphysik II and Center for Collective Quantum Phenomena in LISA+, 
Universitat Tubingen, Auf der Morgenstelle 14, 72076 Tubingen, Germany 

PACS numbers: 74.25.Dw, 74.25.Op, 74.78. Na, 74.78.-w 

Abstract. We investigate the nucleation of superconductivity in a superconducting Al strip under the influence 
of the magnetic field generated by a current-carrying Nb wire, perpendicularly oriented and located underneath 
the strip. The inhomogeneous magnetic field, induced by the Nb wire, produces a spatial modulation of the 
critical temperature T c , leading to a controllable localization of the superconducting order parameter (OP) wave 
function. We demonstrate that close to the phase boundary T c (B ex t) the localized OP solution can be displaced 
reversibly by either applying an external perpendicular magnetic field B ex t or by changing the amplitude of the 
inhomogeneous field. 
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The sensitivity of superconductivity to the local strength of a magnetic field has been exploited during the 
last years to confine superconductivity by applying a non- uniform magnetic field b(r) [TJ [21 E] ■ The experimental 
realization of this "magnetic" confinement can be achieved, e.g., in hybrid superconductor (S) - ferromagnet (F) 
structures and ferromagnetic superconductors. The properties of the ferromagnetic superconductors and the 
S/F hybrids with rather strong exchange interaction between superconducting and ferromagnetic subsystems 
were discussed in the reviews [H O [SJ [Z]- Hereafter we will focus on the flux-coupled hybrids, where the 
interaction between superconducting element and the sources of the magnetic field (e.g., domain walls in the 
ferromagnetic film) occurs via slowly decaying stray fields only j8j (U [10] . 

In general, for thin-film superconducting samples, infinite in the lateral directions, superconducting order 
parameter (OP) wave function first nucleates near the \B z (x, y) \ minima, where B z (x, y) = B ext + b z (x, y) is the 
out-of-plane component of the total magnetic field, B ext is the applied external magnetic field (see arguments, 
e.g., in [10]). Depending on B extl favorable conditions for the appearance of superconductivity can be fulfilled 
either above domain walls in a thick ferromagnetic substrate (domain-wall superconductivity [TTJ [T2J [T3| ) , or 
above magnetic domains of opposite polarity with respect to the B ex t sign (reverse-domain superconductivity 
[TJ1 H21 Uni [TTJ [T5] ) . The external- field-induced crossover between domain-wall superconductivity and reverse- 
domain superconductivity as \B ext | increases can result in an unusual dependence of the superconducting critical 
temperature T c on B ext , which can be nonlinear or even non-monotonous [HI [TjJl [20j [2TJ in contrast to a plain 
superconducting film in a uniform magnetic field. In addition, these domain patterns are periodic in space 
and therefore superconductivity has to be located at all magnetically compensated areas. In order to have a 
singly connected OP solution a non-periodic field profile is needed [22 . It is interesting to note that the stray 
field of a single domain wall in a thin ferromagnetic layer cannot provide domain-wall superconductivity and 
non-monotonous (or, in the other words, reentrant) phase boundary T c (B ext ) due to vanishing of the field at 
large distances from the domain wall [2] . Even if the amplitude of the stray field produced by domain structure 
becomes insufficient to localize superconductivity, such parallel magnetic domains can induce preferential vortex 
motion and giant anisotropy of the critical currents in superconducting films and crystals [23j [24j [25[ [26[ [27j [28] . 




Figure 1. (Color online) Schematic representation of the sample: the Al strip (top gray element) placed 
perpendicular on top of the Nb wire (bottom red element). 

We would like to note that the amplitude of the stray magnetic field and its profile in real S/F bilayers 
are dictated by the saturated magnetization of the ferromagnet and by the period of the domain structure (or 
magnetic dot array), therefore the flexibility of S/F hybrids is limited [HI [20, 29, 30, 3TJ. Full control over the 
amplitude of the inhomogeneous magnetic field can be reached for the hybrid structures, where the ferromagnetic 
subsystem is replaced by current-carrying coils/wires. Such superconductor-electromagnet (S/Em) hybrid 
systems (also called cryotrons) were invented in 1950's and originally considered as superconducting computer 
elements and circuits controlled by local magnetic field of the coils/wires (see classical textbooks and reviews 
[521 l3"3l IM1 133]). However the most transport measurements done on cryotrons were carried out at low 



Localization of superconductivity in superconductor- electromagnet hybrids 



3 



temperatures, when the ability to manipulate the intensive superconducting currents seems to be the most 
effective. To the best of our knowledge the report of the S/Em properties at high temperatures was carried out 
by Pannetier et al. [36] . In this work for the S/Em system, consisting of a plain Al film and a lithographically 
defined array of parallel metallic lines, it was experimentally demonstrated that (i) the T c {B ext ) dependence 
can be non-monotonous for considerably large driving current Io in these metallic lines and (ii) the shape of the 
T c {B e xt) can be reversibly changed as Io varies. 

In this paper we study the influence of the non-periodic magnetic field b^, generated by a single current- 
carrying wire, on the nucleation of superconductivity. The simplicity of this system allows us to fully understand 
the combined influence of the external homogeneous field and the inhomogeneous field by the electromagnet in 
the migration of the superconducting order parameter. To avoid heating effects, this current-carrying wire was 
fabricated from a superconducting material (Nb) with a considerably higher critical temperature T c than the 
investigated microbridge (Al strip). Due to the design of the sample, the perpendicular z— component of the 
magnetic field, playing an important role for thin-film structures, is uniform across the Al strip {y— axis) and it 
varies only along the strip (x— axis), vanishing slowly as the distance from the wire increases. This configuration 
allows us to directly detect the localization of the OP wave function in experiment as the perpendicularly oriented 
external field B ext or the magnitude of the non-uniform field Bq = max \b WjZ (x, y)\ are changed. We also show 
that for low magnetic fields, although the nucleation of superconductivity occurs at those positions, where the 
z— component of the total magnetic field is close to zero, there is still a systematic decrease in T c as a function 
of B ext , resulting from the increasing gradient of the magnetic field at this position of the OP localization. 
This work expands our previous investigation of the vortex dynamics in a similar system [37 by describing the 
influence of the non-uniform field of the wire on the phase boundary T c (B ext ). 

The hybrid samples consist of a 4 /jm wide and 100 nm thick Al strip, patterned by electron-beam 
lithography and lift-off technique, placed perpendicularly on top of a 1.5 /im wide and 50 nm thick Nb wire 
processed by e-beam lithography and Ar ion milling (Fig.[TJ). In order to avoid electrical contact, the Al strip and 
the current-carrying Nb wire are separated by a 120 nm thick insulating Ge layer. All details of the fabrication 
processes are given in Ref. [37 . To investigate the spatial localization of superconductivity in the Al strip, two 
sets of voltage contacts were prepared at distances of 10 /jm and 50 /im from the Nb wire (the inner and outer 
contacts, respectively). 

The normal (N) - superconductor (S) phase boundaries of the Al strip, measured in the perpendicular 
external magnetic field B ext using the outer voltage contacts, are shown in Fig. [U(a) for different currents in 
the Nb wire, I w , while sending a bias current density of 7.5 x 10 3 A/cm 2 through the Al strip. To determine 
the superconducting transition temperature T c a 99% criterion of the normal state resistance R n was used. 

When the control current is zero (I w = 0) and thus the magnetic field is uniform, the resulting phase 
boundary shows the expected linear dependence of T c on B ext [black circles in Fig. [2fa)]. Due to a rather high 
surface to volume ratio of our mesoscopic sample and the high criterion for the determination of T c , it is natural 
to attribute the phase boundary to the appearance of surface superconductivity [38 : 




Applying this equation to the phase transition line T c (B ext ), measured for I w = 0, we determine the 
superconducting critical temperature in zero field T c q ~ 1.26 K, as well as extrapolated to T = the coherence 
length Co ^ 175 nm and the upper critical field B$ = $ Q /[27r£§] ~ 10.4 mT, $ = 2.07 x 10~ 15 Wb is the 
magnetic flux quantum. 

Interestingly, when the non- uniform component of the magnetic field becomes nonzero (I w ^ 0), the 
phase boundaries exhibit a clear enhancement of the critical field for the whole temperature range as shown 
in Fig. [2ja). This enhancement becomes more pronounced as I w increases. Such behavior was observed for 
various S/F hybrids (see, e.g., review [10) and references therein). The described "magnetic bias" is commonly 
explained in terms of the local compensation of the applied magnetic field due to nonuniform component of 
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Figure 2. (Color online) (a) Phase boundaries T c (B ex t), composed using a 99% criterion of the normal state 
resistance of the Al bridge, for different current values in the Nb wire; the zoomed parts of these curves within 
the dashed box are shown in (b). The calculated phase boundaries T c (B ex t) for the same current values I w = 0, 
5, 7.5 and 10 mA are shown in (c). For I w = we used Eq. (1), while Eq. (0 was used for nonzero I w values. 



the field (B"" n — B ext — B if B ext > B ) and trapping of the OP wave function at the locations near the B z 
minima. 

To illustrate the evolution of the superconducting properties in the considered system upon varying T, B ext 
and I w , we performed numerical simulations within the two-dimensional (2D) time-dependent Ginzburg-Landau 
(TDGL) model [35]. For simplicity we assume that the effect of the superfluid currents j s on the magnetic field 
distribution is negligible and consider the internal magnetic field B equal to the field of the external sources 
'Bext + h w (x). This assumptions seems to be valid (i) for mesoscopic thin- film superconductors with lateral 
dimensions smaller than the effective magnetic penetration depth A = A 2 /d (A is the London penetration depth, 
d is the film thickness); (ii) for superconductors at large H and/or T (i.e. close to the phase transition line), 
when the superfluid density tends to zero. In particular, the TDGL equations take the form 

J^+i^V-T (4>~ MV) + (V + iA) 2 ^, (2) 

VV = divj s , j s = - % -t {^* (V + *A) V - c.c.}, (3) 

where ip is the normalized OP wave function, cp is the dimensionless electrical potential |40j . A is the vector 
potential [rot A = (B ext + b WyZ ) z ], T is temperature, r = (1 — T/T c q), u is the rate of the OP relaxation, 
c.c. stands for complex conjugate. For the interfaces superconductor/vacuum or superconductor/insulator the 
boundary condition for ip and tp has the standard form 

s+"-) r *-* (£) r -°- <*> 

where n is the normal vector to the sample's boundary T. 

The self-consistent TDGL modelling for the superconducting sample with realistic dimensions (close to the 
experimental ones) is impossible because of enormous data flow. Therefore, we have restricted our consideration 
and analysis to a mesoscopic superconducting rectangle: length L = 140 an d width W — 20 £o (e.g., if 
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£o ~ 0.175 /im for thin-film Al superconductors, then L ~ 24.5 /im and W ~ 3.5 /im). Since the superconducting 
OP wave function always has maxima near the corners of the sample [HJ |3H S3] > we formally have to assign 
i/j = at x = ±L/2 (both "left" and "right" edges). This simple technical trick guarantees that the finite length 
of the sample and 90°-corners do not strongly affect those OP solutions, which are localized near the control 
wire and thus of main interest in the current study. 

It should be noted that the field induced by the wire is antisymmetric with respect to its middle line x = 0, 
°w,z{x,y) — —b WtZ (—x,y), therefore (i) the effective field compensation favorable for the OP nucleation occurs 
at positive x values (i.e., on the right side of the wire) for B ext > and vice versa; (ii) T c (B ext ) — T c (—B ext ). 
It allows us to consider only positive B ext values without loss of generality. 
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Figure 3. (Color online) Stationary OP wave functions |i/>(x,y)| forming in the superconducting rectangle 
(length 140£o and width 20£o) m the presence of a current-carrying wire (current I w = 10mA, width 9£o)- 
Panel (a) corresponds to T = 1.20 K; panels (b) and (c) are for T = 0.90 K and T = 0.50 K, correspondingly; 
the B ex t values are indicated in the plots. We show only the area to the right of the wire (80 X 20 §o in 
size) in order to better visualize the localized superconducting states trapped at x > 0. The dash vertical lines 
correspond to the projections of the edges of the control wire. White solid lines show the resulting profile of the 
magnetic field B z (x) = B ex t + b w ^ z {x). The condition rjj = at the left and right edges was used for the TDGL 
modelling, Eqs. (2)-(3). 



Figure [3] shows typical stationary OP distributions obtained from the TDGL model for high, intermediate 
and low temperatures (1.20 K, 0.90 K and 0.50 K). 

For finite I w and small B ext the OP wave function expectedly nucleates far away from the control wire 
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Figure 4. (Color online) Phase boundaries calculated for the superconducting rectangle (length 140 £o and 
width 20 £o) in the presence of the current-carrying wire (current I w = 10 mA, width 9 £o) within 2D TDGL 
model [Eqs. (2)-(3)], and ID linearized GL model [Eq. (5)]. We assume the identical parameters for both 
models (lateral dimensions, boundary conditions, vector potential distribution, grid size etc). The two dashed 
lines represent the phase transition line for the development of bulk superconductivity T c 2 = T c o (1 — |/3|) and 
for the critical field of surface superconductivity T c 3 = T c q (1 — 0.59 |/3|), where = B ex t / '-fcS ■ ^ ne svm bols 
show the corresponding points for the images in Fig. \3\ 



(see panel /3 = 0.03 in Fig. G^a), (3 = B ext /B^>), i.e in the area where the field induced by the wire vanishes 
and T c tends to T c0 as |B ea; t| — > 0. At the same moment the OP nucleation near the middle line of the wire, 
where also B z (x,y) « 0, is less energetically favorable due to a larger field gradient dB z /dx at x = 0. The 
depletion of bulk superconductivity at very large distances from the wire by the applied field B ext corresponds 
to a threshold value /3o = 1 — T/T c q ~ 0.048 for T = 1.20 K. The local suppression of the superconducting 
condensate (blue spots) at rather large distances from the wire for (3 — 0.07 indicates the formation of vortices 
by the combined magnetic field of the control Nb wire and the applied magnetic field. For larger B ext the OP 
wave function becomes more localized along the x— axis and trapped near the region where B ext + b WiZ (x) = 0, 
remaining more or less uniform across the Al strip (panels (3 = 0.12 and (3 = 0.18). Thus, for rather high T 
(close to T c0 ) and low B ext , the localized OP solution moves towards the Nb wire as \B ext \ gradually increases, 
until it finally reaches the edge of the current-carrying wire. This process is accompanied by a monotonous 
decrease in T c (see Figs. [5]and|3]), which is a direct consequence of a shrinkage of the typical width of the OP 
solution in the £— direction as |-B ex {| increases (an analog of the quantum-size effect for the Cooper pairs in the 
nonuniform magnetic field |10j). 

For lower T, a completely different evolution of the superconducting properties is observed, when the 
inhomogeneity of the OP wave functions across the strip becomes crucial. Indeed, the formation of the 
localized superconducting state near the minimum of the total field, B™ m « B ext — Bo, leads to the following 
asymptotic behavior at \B ext \ ^ Bq: T* 2 c± T c q x (1 — \B ext — BqI/B^). In addition, a survival of the 
surface superconductivity at large distances, where B z ~ B ext , can be estimated according Eq. (1). Due to the 
different slopes dT* 2 /dB ext and dT c ^/ dB ext and the different offsets, we get the point B* xt « 2.44 Bo, where 
both critical temperatures are equal: T* ~ T c q (1 — 1.44 Bq/B^) (the similar argumentation was presented, 
e.g., in [44]). Substituting B — 0.36 -B^ (estimated for I w = 10mA), we get the threshold temperature 
T* ~ 0.6 K. It means that for intermediate temperatures (T = 0.90 K, panel (b) in Fig. [3]) the gradual increase 
in B ext causes subsequently the development of bulk superconducting state (the plot labelled (3 — 0.20), the 
complete suppression of bulk superconductivity (/3 = 0.40), and then the suppression of surface (or edge- 
assisted) superconductivity at large distances from the wire and the survival of the state localized near the 
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right edge of the wire (f3 = 0.51), since T* 2 > T c3 . Finally, superconductivity survives in a form of 2D patterns 
localized in both directions and centered at the points where the wire's right edge intersects the superconducting 
sample (/3 — 0.60). In contrast to that, for low temperatures (at T = 0.50 K, panel (c) in Fig. [3]) the enhanced 
superconductivity along the wire's edge is suppressed before the destruction of edge-assisted superconductivity 
(images (3 = 0.80 and (3 = 0.90), since T* 2 < T c ^ in this temperature range. 

To show better the observed difference in the shape of localized OP patterns at high and low temperatures 
we calculate numerically the phase transition lines T c (B ext ) using the described 2D TDGL model, Eqs. (2)-(4) 
and ID linearized Ginzburg-Landau (GL) equation 



assuming ip — ip(x,y) and ip — ip(x), respectively, for the same patterns B z and A and the same boundary 
condition ^ = at x = ±L/2. The used ID model was described in detail in [19]. It is easy to see that for low 
Bext and high T both models give almost identical results, what supports our conclusion that the appearing OP 
solutions are almost uniform across the strip and the OP inhomogeneity in the y— direction can be disregarded. 
For intermediate and low temperatures and for large B ext the phase boundaries T c (B ext ) reveal a linear behavior 
with different slopes: dT c /dB ext ~ dT c2 /dB ext for ID model and dT c /dB ext ~ 0.59 dT c2 /dB ext for 2D model. 
It indicates that in the limit \B ext \ ^> B the superconductivity is trapped both near the strip edges (in order 
to correspond to the slope typical for the surface superconductivity) and at the minimum of the local magnetic 
field (in order to explain the parallel shift of the high- field asymptote of T c (B ext ) in higher field). 

Since we are mainly interested in the migration of the OP along the Al strip in low magnetic fields, we can 
propose a very simple description based on ID linearized GL model, neglecting the finiteness of the Al strip 
in the y— direction [2J. If superconductivity is confined within an area where the local magnetic field can be 
approximated by a linear dependence B z (x) = b' w z (xq) (x — xq), where xq is the point of zero total magnetic 
field (B ext + b w . z (xo) = 0), for the phase boundary we obtain as a rough estimate 



Considering the generic case of a non-uniform field b w z — [IqI w / '(2ttx), induced by an infinitely thin cylindrical 
wire carrying a control current I w , at the large distances from the wire one gets \b' w z (xq)\ = 2irBl xt / (p I w ) 
and hence 



which seems to be valid only for low fields, \B ext \ <C Bq, when the OP wave function is located far from the 
current-carrying wire. We compare the experimental data, obtained for low B ext [Fig. HJb)], with the results 
of this model [Fig. HJc)], where we use the same parameters T c o, £o and B^. We conclude that the ID model 
works quite well for describing the T c suppression and the OP migration along the Al strip upon sweeping B ext . 
By changing both B ext and /„, we directly verify that the nucleation of superconductivity is controlled not only 
by the local magnetic field, but also by the gradient of the field in the area where superconductivity is confined. 

It is worth noting that this magnetic field profile is very similar to the field produced by a single domain wall 
in a thin ferromagnetic layer with out-of-plane magnetization [2J. The N-S phase boundary of a superconducting 
film on top of such a ferromagnetic system calculated in [2] looks very similar to that observed in our experiments. 
As a result, we can claim that our S/Em hybrid system behaves as a ferromagnetic film with a pinned straight 
domain wall with tunable saturated magnetization underneath a superconducting thin film, i.e. a situation, 
which cannot be easily achieved in the S/F bilayers. 

Up till now we presented the superconducting properties measured with the outer voltage contacts located 
at the distance 50 /im from the Nb wire. To prove the concept of the OP localization directly, we prepared 
a second pair of inner voltage contacts closer to the Nb wire (10 /an from the wire). The results of these 



tiJj + (V + iA) 2 iJ; = 0, 



(5) 
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Figure 5. (Color online) N/S phase boundaries determined using a 95% criterion of the normal state resistance 
of the Al film for different current values in the Nb wire, measured at the inner contacts of the bridge. 




Figure 6. (Color online) Typical dependences V on B ex t measured for both outer and inner voltage contacts 
at temperatures close to T c o, when one-dimensional OP migration along the Al wire is anticipated. It should be 
noted that normal state resistances R^ n \ = V^ n ) /In — 0.45 Q and = /In ~ 2.53 C, measured 

outer outer/ u inner inner' u ' 

at the corresponding contacts at high magnetic field, differ approximately 5.6 times due to the difference in 
length between these contacts (In = 50 fiA is the bias current). Therefore, for a better visualization we show 
the difference Vouter — Vinner and Vinner + SV, which can be attributed exclusively to the voltage drop outside 
the inner contacts and between the inner contacts, respectively. The offset SV ~ 82 fiV was introduced in order 
to equalize the position of the plateaus on V ou ter — Vinner and Vinner + SV at large B ex t- 

measurements are shown in Fig. [5j In this case to determine the phase boundary T c (B ext ) we used a 95% 
criterion, so that the same condition for detection of superconductivity as the 99% criterion for the 100 /iin 
voltage pad separation is obtained for the inner contacts of distance 20 /im. 

The key finding is that T c measured by the inner contacts is always lower than T c for the outer contacts, 
provided \B ext \ is rather small (up to 0.5 mT). Nevertheless the high-field asymptotes T c (B ext ) for both types 
of contact arrangements expectedly coincides for high fields (compare Fig. EJa) and Fig. [5]). This implies that 
the inhomogeneous superconductivity is located within the inner contacts for high fields and outside these 
contacts for low fields. Indeed, in our magnetoresistive measurements using the inner contacts at rather high 
temperatures and B ext ~ we cannot detect superconductivity nucleating outside these contacts. Therefore the 
measured T c must be lower than the critical temperature for the OP solution localized far away from the Nb wire. 
This is convincing experimental evidence for the field-dependent OP localization in the non-uniform magnetic 
field. We observed that upon increasing B ext , the localized superconductivity shifts toward the wire, leading to 
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a non-monotonous variation of the resistance measured between the inner contacts and, correspondingly, to the 
non-monotonous variation in T c (see Fig. [5]). Unlike the reentrant dependence Vi nner vs. B ext , the voltage drop 
measured at the outer contacts monotonously increases as B ext increases, since localized superconductivity 
always nucleates between the outer contacts and never leaves this area. Considering the inset in Fig. [SI we 
propose the following interpretation of our findings: in the region I the voltage drop V ou t er — Vinner, attributed 
to the area between the outer and inner contacts solely, is minimal and almost independent on B ext , while 
Vinner is maximal, therefore there is no global superconductivity at this temperature in the sample and the OP 
wave function has to be localized between the outer and inner contacts. In the region II both V ou t er — Vi nner 
and Vi nner strongly depend on B ext , and as a result, the OP wave function should be located somewhere in 
the vicinity of the inner contacts. Finally, in the region III the difference Vouter ~ Vinner becomes equal to 
the field-independent value corresponding to the normal state and Vinner is still smaller than its normal value, 
therefore the OP wave function is definitely trapped between the inner contacts (i.e., near the control wire). 
This behavior is in good agreement with the theoretical predictions. 

The position of the T c maximum for the inner contacts measurements can be attributed to the effective 
compensation of the built-in magnetic field at the inner voltage contacts. Due to the small asymmetry in the 
position of the voltage contacts at the opposite sides of the wire as a result of fabrication imperfections, both 
peaks occur at different fields: B ext = — 0.32mT and B ext = 0.26 mT, which is on the order of what we expect 
from the induced field at the voltage contacts for a current of 10 mA. More interestingly, the T c 's corresponding 
to these peaks are slightly different, clearly showing the influence of the field gradient at the point of localization 
on T c . 

Summing up, we have studied the OP localization in an Al strip subjected to an inhomogeneous field with 
tunable intensity induced by a current-carrying wire. The OP migration along the strip upon varying B ext and 
I w has been detected by using multiple voltage contacts. We have shown that the critical temperature at the 
compensated positions is dependent on the local variation of the magnetic field. Interestingly, we demonstrate 
that both reentrant and non-reentrant superconducting phase boundaries can be obtained depending on where 
the voltage drop is recorded. 
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